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Preface 


The  purpose  of  this  study  is  to  use  a  general  technique 
to  analyze  the  radiation  and  scattering  performance  of  a 
logarithmic  spiral  antenna.  This  type  of  antenna  is 
typically  used  in  radar  warning  receivers  because  has 
more  than  one  mode  of  radiation,  which  gives  the  a'cenna  a 
direction  finding  capability. 

The  study  uses  a  relatively  new  general  proyj  known 
as  the  Electromagnetic  Surface  Patch  (ESP)  code  *•-,  analyze 
the  antenna.  Three  antenna  modelling  techniques  are 
explored  and  results  are  compared  to  show  how  tc  .'l  tain  the 
most  accurate  modelling. 

The  first  chapter  gives  a  general  background  of  antenna 
scattering  and  soma  of  the  basic  knowledge  of  spiral 
antennas.  The  second  chapter  gives  a  development  of  the 
method  of  moments,  which  is  a  mathematical  process  used  in 
the  ESP  code.  This  chapter  is  intended  to  give  only  the 
basics  necessary  to  understand  the  moment  method.  The  third 
chapter  introduces  the  three  computer  modelling  schemes. 
The  computer  results  for  the  models  are  given  in  the  fourth 
chapter.  This  chapter  also  compares  the  results  to  draw 
some  specific  conclusions.  The  fifth  chapter  gives  the 
general  conclusions  for  the  study  and  discusses 
possibilities  to  extend  the  study.  Finally,  the  three 
appendices  are  given  to  support  the  computer  modelling 
discussion  in  the  third  chapter. 
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Abstract 

This  study  uses  the  method  of  moments  with  the 
Electromagnetic  Surface  Patch  Code  to  analyze  the  radiation 
patterns  and  currents  developed  on  a  planar  logarithmic 
spiral  antenna.  The  study  also  outlines  a  method  to  analyze 
the  radar  cross  section  of  the  antenna.  Three  computer 
modelling  concepts  are  explored.  The  first  and  most  simple 
is  the  thin  wire  model,  which  is  a  set  of  straight  thxn  wire 
segment?  which  approximates  the  shape  of  the  antenna.  The 
second  concept  e/.ylored  is  that  of  wire  grid  modelling.  The 
purpose  of  thie  model  is  to  check  the  transverse  components 
of  current  along  the  strips  of  the  spiral's  arms.  The  last 
model  is  a  patch  model,  which  is  introduced  to  best  model 
the  structural  scattering  component  of  radar  cross  section. 
Unfortunately,  this  model  proved  to  be  infeasible  unless  a 
very  fast  computer  system  is  used.  Results  are  given  for 
the  eum  and  difference  modes  of  radiation.  These  results 
show  that  the  thin  wire  model  is  adequate  for  radiation 
analysis,  but  the  tips  of  the  arms  must  be  modelled  with  an 
impedance  to  absorb  the  outward  travelling  energy  in  the 
current.  Without  this  loading,  the  arms  will  develop  a 
standing  wave  pattern  for  the  current  distribution,  which  is 
inconsistent  with  experimentation. 


CHAPTER  1 


Introduction 


Background 


One  of  the  major  contributors  to  Lne  radar  cross 
section  of  an  air  vehicle  is  antenna  scattering.  It  is  an 
inherent  and  desired  property  of  an  antenna  to  receive 
and/or  radiate  electromagnetic  energy  efficiently.  Unfortu¬ 
nately,  an  antenna  may  radiate  *>hen  excited  by  incident 
radiation,  as  well  as  when  excited  at  the  feed.  It  is  this 
reradiation  of  energy  which  is  known  as  antenna  scattering 
and  produces  the  radar  cross  section  of  the  antenna. 
Several  antennas  are  exposed  to  incident  radiation  on  any 
conventional  military  airplane.  At  least  one  of  these 
belongs  to  the  radar  warning  receiver,  which  is  designed  to 
intercept  and  identify  incident  electromagnetic  energy.  It 
is  not  a  trivial  task  to  camouflage  this  antenna,  while 
retaining  its  primary  mission  capabilities. 

Conventional  radar  warning  receivers  use  cavity  backed 
flat  spiral  antennas.  These  antennas  are  broadband  and 
polarization  independent  (for  any  linear  polarisation). 
However,  these  antennas  may  have  an  unacceptable  radar  cross 
section . 


Problem 

The  scattering  mechanisms  of  a  radar  warning  receiver 
antenna  must  be  investigated  and  understood  to  reduce  the 


a 


general  method 


radar  cross  section  of  the  antenna.  Also, 
of  analysis  needs  to  be  used  to  determine  the  radiation 
characteristics  of  the  antenna. 

General  Approach 

The  scattering  characteristics  o  a  spiral  antenna  can 
be  derived  from  the  currents  induced  on  the  antenna  by 
incident  electromagnetic  energy.  In  this  study,  the 
reaction  integral  equation  is  used  to  describe  the  interac¬ 
tion  of  the  incident  and  scattered  fields  with  the  current 
distribution  induced  on  the  antenna.  The  reaction  integral 
equation  can  also  be  used  to  describe  the  interaction  of  the 
current  distribution  with  fields  coming  from  other  currents 
that  together  radiate.  A  solution  to  the  reaction  integral 
equation  is  obtained  via  the  method  of  moments. 


Theory 


The 

scattering 

charac 

ter istics 

of  an 

object  are 

normally 

described  by  the  radar  cross 

section 

(RCS)  of  the 

object . 

The  RCS 

is  a 

parameter 

which 

is  directly 

proportional  to  the 

power  returned  to 

a  radar 

which  illumi- 

nates  the 

object . 

The  RCS 
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|es|  =  amplitude  of  scattered  electric  field  at 
observation  point 

jfi!  »  amplitude  of  incident  electric  field  on  object. 
Although  the  RCS  has  units  of  length  squared,  it  is  not 
necessarily  related  to  the  physical  size  of  the  object. 

R.  B.  Green  presented  the  fundamental  concepts  of 
antenna  scattering.  He  started  by  giving  the  simple  circuit 
models  for  a  transmitting  antenna  shovn  in  Figure  1. 


Figure  1.  Equivalent  Circuits  for  Transmitting  Antennas 
Source : [ 2 ] 


The  driving  sources  are  the  current  and  voltage  sources. 
The  nodes  represent  the  terminals  of  the  antenna.  The 
terminal  (input)  current  to  the  antenna  is  denoted  ir  in 
each  case.  The  subscript  r  denotes  radiation.  The  input 
impedance  of  the  antenna  is  Za .  Vt(r?  is  the  voltage 
induced  across  the  terminals  of  the  antenna. 

When  the  terminals  of  a  receiving  antenna  are  short- 


circuited,  the  scattering  antenna  may  be  modelled  with  the 
simple  circuit  in  Figure  2. 


Figure  2.  Short-Circuited  Scattering  Antenna 
Source : [ 2 ] 


For  the  short-circuited  antenna,  the  terminal  voltage,  Vt , 
must  be  zero.  The  short-circuited  terminal  current  is 
denoted  1(0).  The  voltage,  Voc,  is  the  terminal  voltage 
that  would  exist  if  the  terminals  were  open-circuited.  This 
is  easily  seen  because  no  current  would  exist  in  the  circuit 
model  when  the  terminals  are  left  open. 


Mew  a  general  load  impedance  may  be  added  to  the 
antenna  terminals  to  represent  the  input  impedance  of  the 
radar  warning  receiver  circuitry.  This  is  shown  in  Figure 


3.  which  follows. 


J 

Figure  3.  Loaded  Scattering  Antenna  Model 
Source:  [2] 

From  simple  circuit  theory, 

I(Zl)  -  -Voc/(Za  +  Zl)  ( 

From  Figure  2,  it  is  easily  seen  that 

Voc  «  -Za  *  1(0)  ( 

Substituting  equation  (1-3)  into  (1-2), 

I ( Zl )  -  Za  *  I ( 0 ) /( Za  +  Zl)  ( 

Now  the  following  linear  equations  may  be  written  to 
cribe  the  loaded  terminal  voltage  and  current  in  terms 
the  short-circuit  and  radiation  voltages  and  currents. 

I(Zl)  -  Cl  *  1(0)  ♦  C2  *  Ir  ( 


Vt(Zl)  *  Cl  *  vt(0)  +  C2  *  Vt ( r ) 


( 


These  equations  are  justified  because  each  circuit  model  is 
linear  and  any  signal  which  occurs  at  the  terminals  of  the 
loaded  scattering  antenna  (current  or  voltage)  must  be  a 
linear  combination  of  the  terminal  signals  in  the  radiating 


and  short-circuit  models.  Now  Vt(Zl)  =  -Z1*I(Z1),  Vt(r)  » 
Za*Ir,  and  Vt(0)  is  zero.  Therefore,  equation  (1-6)  becomes 

-Z1*I(Z1)  -  C2  *  Za  *  Tr  (1-7) 

Substituting  equation  (1-7)  into  (1-5), 

I(Z1)  -  C1*I(0)  -  Zl* I ( Zl ) /Za  (1-8) 

or  I(Z1)*(1  +  Zl/Za )  »  C1*I(0)  (1-9) 

Substituting  equation  (1-4)  into  (1-9), 

Za*I ( 0  )  * (  1  +  Zl/Za  )/  (  Za  +  Zl)  •»  C1*I(0)  (1-10) 

This  quickly  reduces  to  Cl  ■  1.  Now  from  equation  (1-5), 

Za*I ( 0 ) / ( Za  +  Zl)  -  1(0)  +  C2*lr  (1-11) 

or  (-1  +  Za/{ Za+Zl ) ) *1 ( 0  )  «  C2*Ir  (1-12) 

or  C2  -  -Z1*I ( 0 ) / ( Ir* ( Za  +  Zl))  (1-13) 


A  similar  linear  combination  is  obtained  for  the  scat¬ 
tered  electric  field  from  a  loaded  antenna.  That  is 

Es(Zl)  «  C 1 *Es ( 0 )  +  C2*Er  (1-14) 

where  EsfZl)  *  scattered  electric  field  from  antenna  with 

load,  Zl 

Es(0)  *  scattered  electric  field  from  short- 

circuited  antenna 

fr  ■  radiated  electric  field  from  antenna  with 
a  unit  current  driving  source 

Cl  »  1 


C  2  -  -Z1*I ( 0 )/( Ir* ( Za  +  Zl)) 


Therefore , 

Es(Zl)  =  Es ( 0  )  -  Zl*I(0)*Er/(Ir*(Za+Zl)  )  (1-15) 

Suppose  the  load  is  the  complex  conjugate  of  the  antenna 
impedance.  Then 

fs(Zb)  =  Es ( 0 )  -  Zb*I(0 )*fr/(2*Ra*Ir)  (1-16) 

where  Zb  «  complex  conjugate  of  Za 

Ra  *  real  component  of  antenna  impedance. 

Now  the  scattered  electric  field  from  an  antenna  with  a 
general  load  may  be  written  in  terms  of  the  scattered  field 
from  an  antenna  with  a  conjugate  matched  load.  Substituting 
equation  (1-16)  into  (1-15), 

Es(Zl)  *  Es ( Zb )  +  Zb*  I (0 )*Er/(2*Ra*Ir) 

-  Zl*I(0)*fir/(Ir*(Za  +  Zl)  (1-17) 

which  may  be  rewritten, 

fis(zl)  -  Es ( Zb)  +  (zb-Zl)*I(Zb)*fr/(lr*(Zl+za) )  (1-18) 

where  I(Zb)  ■  Za*I ( 0 ) / ( 2 *Ra ) ,  which  comes  from  equation  (  1- 
4)  when  Zl  ■  Zb.  The  first  term  on  the  right  hand  side  of 
equation  (1-18),  Sa(Zb),  is  known  as  the  structural  mode 
scattering  component  [3:26].  The  remaining  term  on  the 
right  hand  side  of  equation  (1-18)  is  the  antenna  mode  scat¬ 
tering  component  [3:26].  It  completely  vanishes  when  the 
antenna  load  is  conjugate  matched  to  the  antenna  impedance. 

The  term,  (Zb  -  Zl)/(Za  +  Zl),  may  be  defined  as  the 
modified  current  reflection  coefficient,  R  [3:26].  R  is  the 
same  as  a  normal  transmission  line  reflection  coefficient 
when  the  antenna  impedance  is  purely  resistive. 

Suppose  the  structural  mode  scattering  component  for 


some  antenna  is  exactly  zero 


Now  all  scattered  fields 


arise  from  antenna  mode  scattering.  The  definition  in  equa¬ 
tion  (1-1)  can  be  used  to  derive  the  RCS  from  antenna  mode 


scattering. 


lim  4  fr  r 

r-K» 


lim  4  r 
r-K» 


IfiTs  ( zi)l 

M  1 

2  i2  i2 
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(1-19) 


Assume  the  antenna  is  used  as  a  transmitter  and  the  antenna 
load  is  conjugate  matched  to  the  antenna  impedance  (for 
opti:.  um  transmission).  Now  I(Zb)  is  exactly  Ir.  The 
following  manipulations  can  then  be  made  to  define  the 
antenna  mode  RCS  in  terms  of  the  transmitting  antenna  gain, 
G.  The  gain  is  defined  as  (1:37) 

G  »  4*Pi*U/Pin  ( 1 -20  ) 

where  U  -  radiation  intensity 

Pi  -  3. 14159265 
Pin  *  antenna  input  power. 

Since  interest  is  primarily  in  monostatic  RCS  (that  is,  in 

the  direction  of  backscatter ) ,  U  should  be  found  in  the 

direction  of  backscatter.  This  is  simply  (1:33) 

2  .  .2 


U  * 

r  |fr|  /2  *nu 

(  1 

-21  ) 

where 

nu  * 

free  space  impedance. 

The  input  power  to  the  antenna  (for  scattering)  is 

merely 

the  product 

of  the 

incoming  radiation  intensity 

and 

the 

solid  angle 

c  reated 

by  the  effective  aperture 

of 

the 

antenna.  That 

is , 

8 


2  !-  I2  2 

Pin  =  (r  *  j Eij  /2*nu)*(Ae/r  ) 

2  -  I2 

where  (r  *  |e il  /2*nu)  =  incoming  radiation  intensity 


(1-22) 


Ae/r  =  solid  angle  created  by  the  effective 
aperture  of  the  antenna 

Ae  =  effective  aperture  of  the  antenna. 

Substituting  equations  (1-21)  and  (1-22)  into  (1-19), 


o  -  4* n  *  Jr|  *U*Ae/Pin 

Further  substituting  equation  (1-20)  into  (1-23), 


(1-23) 


o  -  IrI  *  G  *  Ae  (1-24 ) 

The  effective  aperture  of  the  antenna  can  be  expressed  as 
[1:61] 

2 

Ae  -  X  *  G/(4  *  tt  )  (1-25) 

where  \  ■  wavelength. 

Therefore , 

2  2  2 

o  -  H  *  X  *  G  /  ( 4  *  it  )  (1-26) 

which  is  the  simplest  expression  for  the  antenna  mode  RCS. 
It  is  interesting  to  note  that  this  expression  completely 
vanishes  when  the  modified  current  reflection  coefficient 
goes  to  zero.  Also,  the  antenna  mode  RCS  is  directly 
proportional  to  the  square  of  the  antenna  gain  in  the 
backscatter  direction.  These  mechanisms  make  sense  because 
an  antenna  will  receive  and  radiate  power  proportional  to 
its  gain  and  reradiation  of  incident  power  is  a  direct 
result  of  the  current  reflection  (for  antenna  mode). 

The  structural  mode  scattering  component  arises  from 
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phi  =  polar  angle  in  the  plane  of  the  spiral 
ro,  a  =  constants. 

Curtis's  approach  was  to  approximate  the  Archimedian  spiral 
as  a  series  of  semicircles.  He  assumed  a  current  distribu¬ 
tion  of  the  sum  of  outward  and  inward  travelling  waves,  with 
a  reflection  coefficient  at  the  open  end  of  the  spiral  to 
account  for  the  inward  travelling  wave.  He  then  proceeded 
to  solve  the  vector  potential  of  a  single  semicircle  in 
closed  form.  Once  this  was  found,  the  radiation  pattern  for 
an  entire  arm  was  given  by  the  sum  of  the  contributions  of 
the  semicircles.  Curtis  made  the  analysis  complete  by 
including  a  second  arm  rotated  180  degrees  physically  from 
the  first  arm  and  driven  with  a  current  at  the  origin  180 
degrees  out  of  phase  from  the  first  arm  current. 

Curtis's  work  is  limited  in  that  it  only  concerns 
Archimedian  spirals  and  assumed  a  current  distribution  on 
the  antenna.  The  current  distribution  he  used  cannot  be 
validated  with  theory  or  experimentation  unless  the  antenna 
is  extremely  small  with  respect  to  a  wavelength. 

Cheo,  Rumsey,  and  Welch  [5]  used  Maxwell's  equations  to 
approach  an  exact  solution  for  the  radiation  pattern  of  an 
antenna  consisting  of  an  infinite  number  of  equally  spaced 
wires  in  the  form  of  coplanar  equiangular  spirals.  This 
type  of  antenna  is  sometimes  known  as  a  logarithmic  spiral 
and  is  characterized  by  the  equation 

-a*phi 

r  =  e  ( 1-28) 

1 1 


where 


r  =  radial  distance  from  origin  to  any  point  on 
the  spiral 

phi  =  polar  angle  in  the  plane  of  the  spiral 
a  *  constant  of  angle  proportionality. 

The  solution  obtained  by  Cheo,  Rumsey,  and  Welch  for  the 
radiated  electric  field  was  a  complicated  integral  expres¬ 
sion  involving  Bessel  functions.  Using  an  asymptotic 

approximation  for  the  Bessel  functions,  they  obtained  a 
somewhat  simplified  expression  for  the  far-zone  electric 
field.  That  expression  is  15: 530] 

a  (1/a) tan  1(a*cos  0) 

|e  i  cos  9  (tan..9/Jl£.;  . — -  (1-29) 

sin  9  \/ 1  +  a^cos20 

for  a  two  arm  logarithmic  (equiangular)  spiral  in  the  x-y 
plane  where 

g  I  -  amplitude  of  phi-component  of  far-zone 

♦I  electric  field 

a  »  constant  of  angle  proportionality 

0  »  angle  from  z-axis 

<{>  *  angle  from  x-axis  projected  on  x-y  plane. 

The  major  limitation  of  their  work  was  that  it  Involved  only 
radiation  of  the  antenna.  Cheo,  Rumsey,  and  Welch  did  not 
predict  the  current  causing  scattering  by  the  antenna. 

K.  K.  Mei  demonstrated  the  feasibility  of  numerically 
solving  an  integral  equation  that  modelled  the  physics  of  a 
radiating  equiangular  spiral  antenna  [6].  He  used  Hallen's 
equation  [7:79]  to  solve  for  the  current  distribution  on 
the  antenna,  which  he  assumed  wa«  composed  of  a  thin  curved 
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wire.  Mei's  technique  for  solving  the  equation  was  essen¬ 
tially  the  method  of  moments,  which  is  discussed  in  detail 
in  Chapter  2.  However,  Mei's  work  may  have  been  in  vain, 
because  Hallen's  equation  assumes  the  current  distribution 
will  have  a  constant  phase  velocity.  This  assumption  proves 
to  be  poor  when  the  wire  has  a  small  radius  of  curvature 
compared  to  a  wavelength  [8].  This  is  to  be  expected 
because  the  equation  was  originally  developed  for  linear 
antennas . 

Joseph  Mosko  discussed  how  a  four  arm  spiral  antenna 
can  be  used  for  direction  finding  applications,  as  in  a 
radar  warning  receiver  [9].  He  showed  that  three  different 
beam-forming  modes  may  be  found  in  the  four  arm  antenna. 
These  modes  were  formed  by  applying  equal  amplitude  voltages 
with  various  phases  to  the  four  feed  points. 

The  primary  mode  is  known  as  the  sum  mode.  It  forms  a 
large  beam  directly  above  the  plane  of  the  spiral.  Xt  is 
excited  by  generating  equal  amplitude  voltages  with  phases 
of  0,  -90,  180,  and  -270  degrees  in  a  clockwise  rotation 
about  the  four  feed  points.  The  sum  mode  creates  a  large 
current  on  the  antenna  around  the  region  where  the  circumfe¬ 
rence  of  a  loop  is  about  one  wavelength. 

The  other  two  modes  are  the  azimuth  and  elevation 
difference  modes.  Together  in  quadrature,  they  form  a  sin¬ 
gle  mode  which  is  simply  known  as  the  difference  mode.  The 
difference  mode  is  excited  by  equal  amplitude  voltages  with 
clockwise  feed-point  phases  of  0,  180,  0,  and  180  degrees. 
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The 

difference  mode 

is  characterized 

by  a 

large  current 

on 

the 

antenna  around 

the  area  where  a 

loop 

circumference 

is 

approximately  two  wavelengths. 

Mosko  stated  that  direction  finding  is  done  by  compa¬ 
ring  the  strengths  of  the  antenna  return  through  the  three 
modes.  This  can  be  done  for  the  passive  antenna  by  using  a 
network  of  phase  shifters  that  may  be  implemented  best  with 
str ipline  [9:91] . 


CHAPTER  2 


THE  METHOD  OF  MOMENTS 

Introduction 

The  method  of  momenta,  or  moment  method,  may  be  simply 
described  as  a  mathematical  technique  used  to  get  an 
approximate  solution  to  an  integral  equation.  Such  an  inte¬ 
gral  equation  may  be  in  the  form  [1:306) 

-  -E1  (f )  (2-1) 

where  f  *  position  vector  in  space 

F' »  position  vector  to  any  point  on  a 
scattering  source. 

This  equation  is  actually  a  form  of  the  'lectric  field  inte¬ 
gral  equation  (EFIE),  which  is  formulated  from  the  boundary 
condition  that  the  total  electric  field  on  the  surface  of  a 
perfect  conductor  (scatterer)  is  zero.  The  term,  E*(F)  ,  is 

the  incident  electric  field  on  the  scatterer.  This  term 
would  usually  be  known.  The  unknown,  I(F'),  i*  the  current 
induced  on  the  surface  of  the  scatterer.  The  term,K(F,F')  , 
is  referred  to  as  the  kernal  of  the  equation,  and  will  take 
on  different  forms  for  different  problem  geometries. 


Expansion 


Functions 


The  method  of  moments  assumes  the  unknown  function  in 


the  integral  equation  may  be  approximated  as  a  linear  combi¬ 
nation  of  some  known  functions.  These  known  functions  are 


called  expansion  functions  or  basis  functions 


They  may 


take  many  forms.  In  general,  the  unknown  current  function 
in  equation  (2-1)  may  be  written  (1:310] 


I(r')  =  1  I„Fn(r')  (2-2  ) 
In  this  form,  n  is  the  subscript  for  the  nth  expansion 
function,  Fn(r')  •  There  are  N  total  expansion  functions. 
As  will  be  seen  later,  allowing  N  to  become  large  may  give  a 
more  accurate  representation  of  the  unknown  current  func¬ 
tion.  i  is  the  nth  expansion  coefficient,  which  may  be  a 
complex  constant. 

There  are  many  forms  an  expansion  function  may  assume. 
The  choice  is  left  to  the  user  of  the  moment  method.  Perhaps 
the  most  simple  is  the  pulse  function.  In  this  case,  the 
prime  coordinate  space  (scattering  surface)  is  broken  into  N 
subsets.  The  nth  expansion  function  can  then  be  described 
as  fo] lows  [1:310]. 


1  for  p'  in  the  nth  subset 


rn(x')  -  ' 

^0  otherwise 


(2-3) 


The  pulse  expansion  function  is  conceptually  simple  to 
understand.  However,  it  is  usually  a  very  inaccurate  repre¬ 
sentation  for  an  unknown  function,  unless  many  subsets  are 
used.  For  example,  suppose  Figure  5  contains  the  actual 
current  distribution  on  a  thin  wire  antenna. 
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Figure  5.  Example  of  Current  Distribution 

As  can  be  seen  in  the  figure,  many  pulse  functions  would  be 
needed  (at  various  heights)  to  accurately  represent  the 
actual  distribution.  To  attempt  to  represent  the  distribu¬ 
tion  with  only  three  pulse  functions,  as  shown  in  Figure  6, 
iu  clearly  a  poor  approximation. 


Figure  6.  Three  Pulse  Approximation  to  Current  Distribution 


For  the  previous  example,  perhaps  as  many  as  one  hundred 


pulse  functions  would  be  desired  for  an  accurate  representa¬ 
tion  of  the  function.  The  number  of  expansion  functions 
needed  to  accurately  represent  the  function  can  be  substan¬ 
tially  reduced  by  the  choice  of  expansion  functions.  For 
example,  a  triangle  function  may  be  defined  such  that 

-  1  for  (r.-l)a-2a  <  x  <  (n-l)a 

Fn(x)  *  0  otherwise.  (2-4) 
For  this  triangle  function,  a  single  expansion  function 
( n*  1  )  with  a.i  expansion  coefficient  of  one  would  exactly 
represent  the  distribution. 


The  Piecewise  Sinusoidal  Expansion  Function 

One  of  the  more  popular  expansion  functions  used  in  the 
method  of  moments  is  the  piecewise  sinusoidal  function.  It 
is  characterized  by  the  following  equation  [1t323}> 


VX> 


(  sin6(x-x^  ,) 
n-l 

sin6(x  -x  . ] 
n  n-l 


sin8(xn+1-x) 

sin8(x  , -x  ) 
n+l  n 


for  x  <  x  <  x 
n-l  n 

tor  x„  <  x  <  xn+1 


0  otherwise 


and  0  «  27I/X 

X  =  wavelength.  (2-5) 

This  function  is  shown  on  the  next  page  in  Figure  7. 

The  piecewise  sinusoidal  expansion  function  is  usually 

spaced  evenly  so  that  x  ,  -  x  is  a  constant  for  all  n. 

c  n+l  n 


Figure  7.  Piecewise  Sinusoidal  Function 
Source :  [  1  :  324] 

Furthermore,  the  functions  overlap  when  they  represent 
unknown  distribution.  This  is  shown  below  in  Figure  8. 


Figure  8.  Overlapping  Piecewise  Sinusoidal  Functions 
Source : [1:324] 


When  x  -  x  is  small  compared  to  a  wavelength 


the  piecewise  sinusoidal  function  approaches  the  triangle 
function  described  earlier. 


resting  Functions 


Once  the  unknown  function  in  an  integral  equation  is 
written  as  a  linear  combination  of  expansion  functions  with 
unknown  coefficients,  a  testing  procedure  can  be  applied  to 
convert  the  integral  equation  into  a  set  of  linear  equa¬ 
tions.  For  example,  consider  the  integral  equation  (2-1) 
with  the  expansion  functions  of  equation  (2-2)  substituted 
for  the  current  distribution.  That  is 


fZl  F  (r' )K(r,r' )dr' 
n  n 


-E1 (?) . 


(2-6) 


A  general  testing  procedure  is  to  multiply  each  side  of 
equation  (2-6)  by  some  testing  (or  weighting)  function  and 
integrate  over  the  unprimed  coordinate  space.  The  summation 
and  double  integration  can  then  be  interchanged  to  give 


Zln//Fn(r’)Gm(r)K(f,r,)dr'dr  =  -/e1 (r)Gm(r)dr 


where  g  (r)  *  the  rath  testing  function, 

m 


(2-7) 


Suppose  N  testing  functions  are  used.  Now  there  are  N 


equations,  each  of  the  form 


Si  z  =  v 

n  mn  m 


(2-8) 


where 


Z  =  //F_(r')G  (r)K(r,r')dr'dr 

mn  11  111 


V  =  -/E1(r)Gm(r)dr. 

m  * 
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In  the  application  of  the  method  of  moments, 


all  terms 


necessary  to  find  Zmn  and  Vm  would  either  be  known  or 
assumed  (the  testing  and  expansion  functions).  The  problem 
is  then  reduced  to  a  matrix  problem.  This  can  be  written  as 

[Z] [I]  -  [V]  .  (2-9) 
[Z]  is  the  impedance  matrix  with  entries,  Zmn,  for  the  mth 
row  and  nth  column.  It  is  a  square  matrix  of  order  N.  [I] 
i 8  a  single  column  matrix  with  the  coefficient,  In,  in  the 
nth  row.  [V]  is  a  single  column  matrix  with  Vm  in  the  mth 
row.  [I]  and  [V]  each  have  N  entries. 

To  solve  the  matrix  problem,  the  impedance  matrix  must 
be  inverted.  The  solution  is  written  notationally  as 

[I]  =  (Z] _1  [V]  (2-10) 
With  the  coefficients.  In  for  n  -  1  to  N,  found  from  equa¬ 
tion  (2-10),  the  unknown  current  distribution  on  the  scat¬ 
tering  surface  is  determined  by  equation  (2-2). 

The  solution  from  the  moment  method  is  not  necessarily 
the  true  current  distribution  on  a  scattering  surface.  It 
is  only  an  approximation  because  the  true  continuous  inte¬ 
gral  equation  representing  the  physics  is  approximated  by  a 
set  of  discrete  linear  equations.  The  approximation  becomes 
better  when  more  linear  equations  are  used  to  replace  the 
integral  equation.  Hence,  to  get  a  better  solution,  more 
expansion  and  testing  functions  must  be  used. 

Another  approach  to  improving  the  moment  method  solu¬ 
tion  is  to  make  intelligent  choices  for  expansion  and  test¬ 
ing  functions.  As  seen  in  Figures  5  and  6,  the  pulse  func- 
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tion  is  not  a  very  intelligent  choice  for  an  expansion  func¬ 
tion  in  some  cases,  unless  many  pulse  functions  are  used. 
Using  many  pulse  functions  may  seem  like  a  reasonable 
approach,  but  this  would  create  a  large  impedance  matrix, 
which  would  be  difficult  to  invert.  In  fact,  the  amount  of 
computer  time  necessary  to  invert  the  impedance  matrix  and 


solve 

equation 

(2-10 ) 

is  proportional  to 

N3 

[1:340]. 

The  computation 

time 

consideration  makes 

an 

intelligent 

choice 

of  expansion  and  testing  functions 

to 

reduce  the 

order  of  the  impedance  matrix,  an  attractive  approach  to  a 
real  moment  method  problem. 

Point  Matching 

One  of  the  easiest  ways  to  implement  the  method  of 
moments  is  to  use  point  matching.  This  is  just  the  name  for 
the  procedure  of  using  delta  functions  for  testing  func¬ 
tions.  These  are  defined  below. 

S(x-a)  =  0  for  all  x  j  a  (2-11) 

/C6(x- a)  dx  =  1  for  b<a<c  (2-12) 

b 

0  otherwise 

The  delta  function  also  has  the  property  that  [10:273] 

£  6  (x-a)  f  (x)dx  =  f(a)  (2-13) 

for  any  function  f(x). 

Point  matching  is  used  when  N  expansion  functions  are 
chosen  so  that  each  expansion  function  is  non-zero  over  one 
of  N  subsets  of  the  prime  coordinate  space  (scattering  sur¬ 


face)  . 


These  subsets  may  overlap,  but  only  one  expansion 


function  can  be  non-zero  at  a  particular  point  on  each 


subset  where  point  matching  will  be  used.  Either  the  pulse 
or  piecewise  sinusoidal  functions  may  be  used  for  the  point 
matching  procedure. 

For  an  example  of  point  matching,  suppose  a  straight 
wire  along  the  x-axis  becomes  a  scatterer  in  a  problem.  The 
wire  may  be  divided  into  N  non-overlapping  subsets  and  pulse 
expansion  functions  may  be  used  for  the  induced  current 
distribution  on  each  subset.  Point  matching  may  be  applied 
to  the  median  location  in  each  subset.  These  points  would 
be  designated  xn,  where  n  *  1  to  N. 

Using  the  electric  field  integral  equation  in  equation 

(2-1),  the  moment  method  gives  the  entries  to  the  impedance 

and  voltage  matrices  in  equation  (2-8).  Now  the  rath  testing 

function,  g  (r)  »  is  just  the  delta  function  at  xn.  That  is 
m 

Gm(r)  ■  6  (x-xm) .  (2-14) 

Using  the  property  of  equation  (2-13),  the  voltage  matrix 
entries  become 


V  =*  -/e*  (x)  5  (x-xm)dx 

m 

■  -E1(xm) 

(2-15) 


and  the  impedance  matrix  entries  become 

Z  =  //f_  (x ' )  6  (x-xm)  K  (x ,  x  ‘ )  dx '  dx 

mn  n 

=  /Fn(x'  )K(xm,x'  )dx'  .  (2-16) 

The  main  advantage  of  point  matching  is  that  entries 
for  the  voltage  and  impedance  matrices  may  be  found  easier 
than  if  any  other  testing  function  were  used.  This  is  a 


direct  consequence  of  the  property  in  equation  (2-13),  which 
removes  one  of  the  integrations  required  for  computing  each 
of  the  entries.  This  advantage  makes  point  matching  an 
attractive  way  to  implement  the  method  of  moments  when  cal¬ 
culations  are  done  by  hand.  However,  if  numerical  integra¬ 
tion  can  be  performed  on  a  computer,  point  matching  does  not 
offer  significant  advantages  over  other  methods. 

Point  matching  ensures  that  the  original  integral  aqua¬ 
tion  representing  the  physics  of  a  problem  is  exactly  satis¬ 
fied  at  the  point  matching  locations.  What  happens  between 
these  points  is  determined  by  the  choice  of  expansion  func¬ 
tions.  As  in  all  moment  method  procedures,  accuracy  in  the 
solution  will  greatly  increase  when  more  expanse  in  functions 
are  used.  In  the  case  of  point  matching,  this  will  mean  the 
points  where  the  solution  is  exact  (integral  equation  is 
enforced)  are  close  together. 

Galerkin ' 8  Method 

Another  popular  procedure  used  to  implement  the  method 
of  moments  is  Galerkin' s  method.  This  method  can  simply  be 
described  as  the  procedure  of  using  the  same  functions  for 
both  expansion  and  testing  functions. 

One  of  the  interesting  applications  of  Galerkln's 
method  is  to  use  piecewise  sinusoidal  expansion  and  testing 
functions.  Experience  has  shown  that  this  procedure  is  both 
numerically  efficient  and  highly  accurate  [1:324].  It  is 
also  possible  to  use  Galerkin’s  method  with  other  expansion 
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functions,  such  as  the  pulse  function.  However,  experience 
has  shown  that  the  additional  accuracy  obtained  by  the 
pulse-pulse  Galerkin  procedure  versus  the  pulse  point  match¬ 
ing  procedure  is  not  worth  the  extra  effort  [1:324]. 

In  general,  Galerkin's  method  will  result  in  the  mini¬ 
mum  mean  squared  error  for  a  given  number  and  type  of  expan¬ 
sion  function  [2]. 

Reaction  Integral  Equation 

The  reaction  integral  equation  is  a  convenient  expres¬ 
sion  relating  the  reaction  of  various  fields  to  each  other 
for  many  problems.  To  understand  the  development  of  the 
equation,  the  concept  of  reaction  must  first  be  understood. 
This  concept  was  developed  by  V.  H.  Rumsey  in  1954. 

Suppose  there  are  two  electromagnetic  sources  in  an 
isolated  universe.  These  sources  can  be  denoted  by  the 
names,  a  and  b.  Suppose  source  a  contains  the  electric  and 
magnetic  currents  of  dJ(a)  and  dK(a),  respectively,  on  some 
differential  volume  in  the  source.  I!  (  a )  and  TJ(a)  can  repre¬ 
sent  the  electric  and  magnetic  fields  generated  by  source  a. 
Similarly,  1(b)  and  H(b)  represent  the  fields  generated  by 
source  b.  Now  the  reaction  of  sources  a  and  b,  <a,b>,  can 
be  defined  below  [11:1483]. 

<  a ,  b>  =  ///[  I(b),dT(a)  -  H  (b) -die  (a)  ]  (2-17) 
This  integration  is  done  over  the  volume  of  source  a. 

Using  the  reciprocity  theorem  [12:477],  it  follows  that 


<  a .  b> 


<b ,  a> 


(2-18) 


as  long  as  the  sources  have  finite  volumes  and  the  medium 
between  the  sources  is  isotropic. 

The  concept  of  reaction  can  now  be  applied  to  a  moment 
method  problem.  Suppose  a  thin  wire  is  a  scatteicr  with  a 
plane  wave  incident  to  induce  an  unknown  current  on  the 

wire.  The  wire  can  be  broken  into  a  finite  number  of  seg¬ 

ments  with  pulse  expansion  functions  representing  the 
unknown  current  distribution.  Let  s  denote  the  source  of 
scattered  fields  and  i  denote  the  source  of  incident  fields. 
Assuming  the  thin  wire  is  a  perfect  conductor,  the  fields 
from  s  and  i  should  equal  zero  on  the  surface  of  the  wire. 
It  is  desired  to  have  no  reaction,  or  coupling,  between  the 

incident  and  scattered  fields  with  a  test  source  (the  mth 

testing  function,  denoted  m).  Therefore 

<s+i,m>  «  0.  (2-19) 

This  can  be  directly  transformed  by  a  property  of  reaction 

[11:1484]  into 

<s,r>  +  <i,m>  ■  0.  (2-20) 

Using  the  definition  of  reaction,  equation  (2-17),  this 

becomes 

-  Hjn'dK3]  +  /// [Em* d^-iydK1]  =  0  (2-21) 

which  is  the  reaction  integral  equation  (RIE). 

If  a  set  of  electric  or  magnetic  test  sources  are  used, 
the  reaction  integral  equation  reduces  to  the  electric  and 
magnetic  field  integral  equations,  respectively  [13:6]. 

To  obtain  a  more  convenient  form  of  the  RIE  for  the 
method  of  moments,  the  property  of  equation  (2-18)  can  be 
applied  to  the  second  term  of  equation  (2-21).  Therefore 


(  2-22  ) 


fill Em'd3S  “  Hm*dKs]  +  Iff tEi-dJm  -  H1’^]  =  0. 

Assuming  a  scatterer  is  almost  a  perfect  conductor,  the 
magnetic  current  causing  scattering  and  the  magnetic  current 
on  a  test  source  (on  the  body)  can  be  neglected.  Now  the 
RIE  reduces  to 

///fm-d5s  +  ///Ei-dJm  =  0.  (2-23) 

The  unknown  in  this  equation  is  the  current  on  the  surface 
of  the  scatterer,  dJs  ,  which  may  be  more  conveniently 
written  as  3s dV  •  Equation  (2-2)  can  be  used  to  repre¬ 
sent  the  unknown  current  distribution.  Therefore, 

Zln//£m-PndS  +  iml‘63n  =  0.  (2-24) 

Note  that  it  has  been  assumed  that  all  scattering  current 
will  be  on  the  surface  of  the  scatterer.  The  voltage  and 
impedance  matrix  elements  may  then  be  defined  as  follows. 

Vm  -  imi-d3m  (2-25) 

Zmn  -  -///£m-FndV  (2-26) 

These  elements  may  be  easily  calculated  once  expansion  and 
testing  functions  are  assumed,  and  the  electric  field  from 
the  rath  testing  function,  Em  ,  is  derived.  This  electric 
field  may  be  computed  by  the  standard  radiation  integral 
approach  [1:25]  on  the  source,  Jm  .  The  solution  to  the  E 
field  for  a  piecewise  sinusoidal  testing  function  on  a  thin 
wire  is  known  in  closed  form  (14:368-370]. 

The  form  of  the  RIE  in  equation  (2-23)  can  still  be 
used  if  the  conductor  is  net  perfect,  as  long  as  it  is  a 
reasonably  good  conductor  so  that  the  magnetic  current 
induced  is  small.  The  only  modification  to  the  previous 
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procedure  is  to  add  a  small  impedance  term  to  the  diagonal 
entries  in  the  impedance  matrix.  This  addition  is  simply 
the  actual  surface  impedance  on  the  mth  subset  (for  the  Zmm 
entry)  of  the  scatterer.  This  is  the  same  impedance  matrix 
modification  used  in  lumped  loading  [1:345]. 

One  of  the  main  advantages  of  the  RIE  over  the  EFIE  in 
equation  (2-1)  is  the  evaluation  of  the  diagonal  impedance 
matrix  entries,  Zmm.  In  the  EFIE,  these  terms  become 
unstable,  because  the  kernal  of  the  equation  (the  free  space 
Green's  function)  has  a  singularity.  This  singularity  is 
avoided  in  computing  Zmm  from  equation  (2-26). 


The 

solution 

to  the  reaction  integral 

equation 

only 

gets  the 

current 

that  causes  scattering. 

Once  this 

is 

found,  the  actual 

scattered  fields  from  the 

surface  may 

be 

obtained 

by  the 

standard  radiation  integral  approach 

of 

antenna  theory  [1:25].  The  radar  cross  section  of  the  scat¬ 
terer  can  then  be  calculated  using  equation  (1-1)  with  the 
range,  r,  some  large  finite  number  (unless  the  electric 
field  i 8  found  by  exact  integration,  in  which  case  the  range 
dependence  will  cancel  .a  the  RCS  formula) . 

Radiation  Problems 

Although  the  entire  chapter  thus  far  has  been 
devoted  to  applying  the  method  of^moments  to  scattering 
problems,  the  radiation  problem  of  an  antenna  may  also  be 
approached.  For  this  problem,  the  moment  method  only  finds 
the  current  distribution. driven  on  the  antenna  for  radia- 
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tion . 


The  only  difference  between  radiation  and  scattering 
problems  is  the  source  of  the  current  distribution.  For 
scattering  problems,  the  source  is  an  incident  electric 
field.  As  shown  in  equations  (2-8 ),( 2-25 )  ,  and  (2-26),  the 
incident  field  is  used  for  voltage  matrix  entries  only. 
This  suggests  that  the  impedance  matrix  entries  are  indepen¬ 
dent  of  the  current  distribution  source,  which  is  true.  The 
same  impedance  matrix  is  used  for  both  radiation  and 
scattering  problems  on  the  same  antenna.  This  makes  sense 
because  the  impedance  matrix  entries  are  only  mutual  impe¬ 
dances  between  different  subsets  of  the  antenna. 

There  are  several  possible  source  models  that  can  be 
used  for  radiation  problems.  The  most  popular  for  wire 
antenna  radiation  is  the  delta  gap  model  [1:336]  which  is 
shown  below  in  Figure  9. 


— >J  6  — 


Figure  9.  Delta  Gap  Model 
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What  is  shown  in  Figure  9  is  one  subset  of  a  wire  antenna  in 
which  a  delta  gap  generator  is  found.  The  delta  gap  genera¬ 
tor  maintains  a  finite  voltage  over  an  extremely  small  gap. 
The  same  finite  voltage  is  found  across  the  entire  wire 
subset.  The  delta  gap  generator  has  the  advantage  of  making 
the  calculation  of  the  voltage  matrix  entries  trivial.  Each 
voltage  matrix  entry  is  simply  the  voltage  across  the 
corresponding  segment  of  wire.  Therefore,  each  entry  is 
zero  except  for  those  for  segments  which  have  delta  gap 
generators,  in  which  case  the  finite  voltage  is  the  entry. 

Regardless  of  how  the  sources  are  modelled  in  a  radia¬ 
tion  problem,  the  moment  method  procedure  is  to  multiply  the 
inverted  impedance  matrix  and  the  source-obtained  voltage 
matrix.  This  will  give  the  coefficients  for  the  expansion 
fu'.  one  representing  the  current  distribution.  This 
current  distribution  can  be  translated  into  the  radiated 
field  pattern  through  the  standard  radiation  integral  tech¬ 


nique  [1:25]. 


CHAPTER  3 


Computer  Modelling 

This  project  makes  use  of  the  Electromagnetic  Surface 
Patch  (ESP)  Code,  Version  II  -  Polygonal  Plates  and  Wires. 
This  code  was  chosen  for  its  versatility  in  creating  a 
quadrilateral  patch  antenna  mnd  ?  i  which  accurately  models 
the  structural  mode  scattering.  This  is  further  explained 
later  in  this  chapter. 

The  antenna  to  be  modelled  in  this  study  is  a  four  arm 
planar  equiangular  (logarithmic)  spiral  composed  of  flat 
curving  strips.  The  strips  flange  outward  as  the  radius  of 
the  spiral  increases  to  help  preserve  the  frequency 
independence  of  the  antenna.  The  feed  to  the  four  arms  of 
the  antenna  is  hidden  inside  an  electromagnetic  shield. 

Thin  Wire  Modelling 

The  first  and  most  simple  model  considered  for  the 
antenna  is  the  thin  wire  model.  This  model  consisted  of 
straight  thin  wire  segments  to  approximate  the  median  curve 
of  the  arm  strips.  Straight  segments  are  used  because  the 
exact  radiation  pattern  from  a  straight  segment  with  a 
piecewise  sinusoidal  distribution  of  c  i.'tunt  is  known.  The 
first  thin  wire  model  used  was  a  two  arm  model,  in  which  two 
opposite  arms  of  the  four  arm  antenna  were  left  out.  This 
is  shown  on  the  next  page  in  Figure  10. 


rxyuxtj  i  i  •  nun  wire  noaexxmg 
Figure  11  also  shows  the  way  the  thin  wire  modelling 
approximates  the  median  curve  of  the  arm  strips.  In  the 
original  modelling  with  thin  wires,  the  ends  of  each  of  the 
arms  was  left  as  an  open  circuit.  That  is,  the  out-most 
segment  of  the  arm  merely  terminated  in  free  space.  Wire 
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This  feed  model  offered  geometrical  symmetry  and  the 
versatility  of  quickly  interchanging  feed  voltages  to 
produce  different  modes  of  radiation.  The  feed  voltages  to 
produce  the  sum  and  difference  modes  of  radiation,  as 
described  in  the  first  chapter,  are  given  below  in  Table  I. 


Table  I 
Feed  Voltages 


Feed  Arm 

Sum  Patt 

Vol  .  (V) 

Diff.  Patt.  Vol 

A 

1 

Ll 

1 

t±° 

B 

1 

/  9  cP 

1 

/l  8  0  ° 

C 

'L 

'  180° 

1 

Zjl° 

D 

'L 

270° 

1 

/  1  8  0  ° 

(  four 


subroutine  listing  the  geometry 
arms)  is  given  in  Appendix  A. 


input  for  this  model 


Wire  Grid  Modelling 

A  wire  grid  model  was  developed  to  check  if  transverse 
currents  could  be  formed  on  the  arms  to  significantly  affect 
the  radiation  patterns.  Figure  13,  on  the  next  page,  shows 


the  concept  of  wire  grid  modelling  of  the  arms 


Ground  Point 


Figure  13.  Wire  grid  Modelling 
As  seen  in  the  figure,  the  wire  grid  model  is  a  set  of 
straight  wire  segments  connecting  points  chosen  on  the  inner 
and  outer  edges  of  the  arm  strips.  The  points  are  chosen  at 
polar  increments  of  45  degress  by  using  emperical 
logarithmic  equations  derived  to  closely  fit  the  actual 
strip  edges.  TK“  polar  increment  angle  (45  degrees)  was 
increased  from  the  thin  wire  model  (in  which  it  was  10 
degrees)  because  the  wire  grid  model  used  more  expansion 
modes  with  its  cross  hatches. 

The  emperical  equations  for  the  inner  and  outer  edges 
of  the  strips  are  given  as  follows. 


It  should  be  noted  that  the  ground  point  is  at  the 
origin  of  a  Cartesian  coordinate  system  while  the  wire  grid 
model  lies  in  the  X-Y  plane.  This  is  similar  to  the 
geometry  of  the  thin  wire  model.  A  subroutine  listing  for 
the  geometry  input  of  this  model  is  given  in  Appendix  B. 

Patch  Modelling 

The  third  model  considered  is  the  patch  model.  It 
consists  of  thin  quadrilateral  coplanar  plates  which  conform 
to  the  shape  of  the  actual  arm  strips.  The  corners  of  the 
quadrilaterals  are  the  same  as  the  points  used  in  the  wire 
grid  model.  The  centers  of  the  segments  of  wire  in  the  wire 
grid  model  are  the  edges  for  the  quadrilateral  patches. 
Therefore,  the  wire  grid  modelling  in  Figure  13  also 
demonstrates  the  geometry  of  patch  modelling. 

The  success  of  patch  modelling  is  highly  dependent  on 
the  ability  to  calculate  the  mutual  impedance  terms  that 
determine  the  impedance  matrix.  The  difficulty  in  finding 
the  patch  to  patch  impedances  constrains  the  computation 
time  dramatically.  This  is  very  different  from  the  cases  of 
thin  segment  impedances,  which  can  be  computed  by  a  line 
integral  of  a  simple  dot  product  (ref.  eq.  2-26).  In  these 
cases  (e.g.  thin  wire  modelling  and  wire  grid  modelling)  the 
impedance  matrix  inversion  (ref.  eq.  2-10)  dominates  the 
computation  time  for  the  computer. 

In  order  to  understand  the  impedance  matrix  entries  of 
the  patch  model,  th  expansion  functions  should  first  be 
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studied 


The  expansion  function  for  patches  used  in  the  ESP 


code  is  a  modified  version  of  the  piecewise  sinusoidal 
function.  It  is  defined  by 


The  t 
which 


,  *  sin  kl 

Js  =  cl 

sin  kLu) 

where  k  *  2v/\ 

erms  in  equation  (3-6) 
gives  the  surface  patch 


are  illustrated 
geometry  . 


in 


(3-6) 


Figure  15, 


R) 


l/L 


Figure  15.  Surface  Patch  Geometry 
Source :  [13:12] 

The  constant,  c,  in  equation  (3-6)  is  chosen  so  that  the 


expansion  function  has  an  amplitude  of  one  ampere  all 
the  edge  of  the  patch  marked  with  a  double  line. 


The  impedance  entry  calculation  involves  the  surface 
integral  of  a  dot  product  of  the  modified  expansion  function 
in  equation  (3-6)  with  the  electric  field  radiated  by 
another  patch  with  a  similar  expansion  function  representing 
its  current  distribution.  There  are  two  reasons  why  this 
impedance  calculation  is  much  more  difficult  than  wire  to 
wire  impedance  calculations.  First,  integration  must  be 
done  over  two  dimensions  instead  of  one.  And  second,  the 
radiated  electric  field  from  this  modified  piecewise 
sinusoid  is  not  as  simple  as  the  radiated  electric  field 
from  a  thin  wire  segment. 

To  illustrate  the  enormous  amount  of  computation 
required  to  solve  patch  model  problems,  the  following 
experiment  was  devised.  The  ESP  code  was  used  to  produce 
radar  cross  section  results  for  some  patch  models.  The 
first  model  was  simply  one  patch  conforming  to  the  first 
(most  inward)  patch  of  one  arm  of  the  antenna  (ref.  Fig. 
1*/.  The  second  model  was  the  first  two  patches  of  one  arm 
of  the  antenna.  The  final  models  used  five  and  ten  patches 
of  one  arm.  Table  II,  on  the  next  page,  gives  the  CPU 
(central  processing  unit)  computation  time  for  a  VAX  11/750 
computer  to  get  similar  RCS  data  on  each  of  these  models. 
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Table  II 


#  Patches 
1 
2 
5 

10 


Patch  Model  Computation  Time 
in  Model  CPU  Time 

1  minute 
6  minutes 
46  minutes 
8  hours 


To  adequately  model  the  entire  antenna,  88  patches  are 
necessary  (at  45  degree  polar  angle  increments).  Since  the 
computation  time  grows  almost  exponentionally  with  the 
number  of  patches,  a  very  conservative  estimate  of  the  time 
required  to  get  RCS  data  of  the  full  antenna  is  over  two 
weeks  of  CPU  time.  With  time  sharing,  the  actual 
calculation  could  be  much  lonqer. 

It  is  a  shame  that  fuli  patch  modelling  cf  the  antenna 
is  not  feasible,  since  it  is  believed  the  patch  model  would 
most  accurately  represent  the  structural  mode  scattering  of 
the  antenna's  RCS. 

An  explaination  of  how  the  main  program  of  the  ESP  code 
is  altered  to  handle  patch  modelling,  along  with  a  sample 
patch  generating  subroutine,  is  given  in  Appendix  C. 


41 


Chapter  4 


Results 

The  first  model  used  for  experimentation  was  the  two 
arm  thin  wire  model  shown  in  Figure  10.  The  purpose  of  this 
experiment  was  to  verify  the  lower  frequency  cutoff  of  1  GHz 
for  the  spiral  antenna.  This  corresponds  roughly  to  the 
frequency  where  one  wavelength  equals  the  outer 
circumference  of  the  antenna.  Table  III  gives  the  maximum 
gain  (along  the  z  axis  as  the  antenna  sits  in  the  X-Y  plane) 
of  the  theta  and  phi  components  of  radiation  of  this  model 
for  various  frequencies.  In  this  study,  the  theta  and  phi 
components  are  from  the  standard  spherical  coordinate  system 
with  theta  being  the  angle  from  the  x  axis  of  a  projection 
on  the  X-Y  plane  [15:364]. 

Table  III 

Gain  of  Components  in  Two  Arm  Model 
Frequency  (MHz)  Theta  Comp.  Gain  (dB)  Phi  Corap.  Gain  (dB) 


250 

1.094 

-7.527 

500 

1.031 

-4.705 

1000 

-2.101 

1 .099 

1500 

-6. 307 

3.  173 

2000 

1.663 

0.944 

3750 

3.043 

2.404 

To  check  the  bandwidth  of  this  circularly  polarized  antenna, 
it  is  better  to  determine  whether  the  antenna  is  polarized 


left  hand  or  right  hand  and  then  find  the  maximum  gain  of 
that  component.  This  can  be  done  by  taking  the  square  root 
of  the  sum  of  the  squares  of  the  components  of  the  fields. 
To  do  this,  the  field  first  must  be  found  by  converting  gain 
from  decibels  by  the  formula 

Gain  -  20  *  In  ( FS  )  (4-1 ) 
where  FS  *  field  strength. 

This  it.  done  to  get  the  field  strength  of  the  components, 
and  once  the  field  strength  of  circular  polarization  is 
found,  the  formula  can  be  used  to  reconvert  into  decibels. 
The  result  is  shown  in  Table  IV. 


Table  IV 


Gain  of  Circular  Polarization  (Two-Arm  Model) 


Frequency  (MHz) 
250 
500 
1000 
1500 
2000 
3750 


Circ.  Pol.  Gain  (dB) 
4.62 
5.50 
6.56 
6.45 
8.24 
9.66 


Assuming  the  maximum  gain  is  about  9.66  dB,  the  antenna 
bandwith  cutoff  gain  should  be  around  6.66  dB,  which  is  3  dB 
lower.  This  is  almost  exactly  what  happens  at  1  GHz. 
Therefore,  the  cutoff  limit  of  a  wavelength  fitting  the 
outer  circumference  has  been  demonstrated. 


Since 


the 


validity  of  the  ESP  code 


has 


been 


demonstrated  with  the  two  arm  model,  the  four  arm  thin  wire 


model  was  used  for  experimentation.  The  sum  patterns  for 
the  phi  and  theta  components  ^f  polarization  are  shown  for  a 
frequency  of  2  GHs.  in  Figures  16  and  17.  The  maximum  gain 
for  the  phi  and  th.te.  polarizations  are  3.328  dB  and  3.330 
dB,  respectively.  Thw  bandwidths  (3  dB)  of  the  phi  and 


Figure  17. 


Figures  16  and  17  seem  to  suggest  that  the  linear 
polarization  independence  of  a  four  arm  spiral  breaks  down 
for  large  angles  off  the  antenna  boresight  (the  z  axis). 

The  difference  patterns  for  the  phi  and  theta 
components  of  radiation  at  2  GHz  is  shown  in  Figures  18  and 


Figures  18  and  19  also  show  a  dramatic  breakdown  in  the 


polarization  independence  beyond  60  degrees  from  the  z  axis. 
The  theta  polarization  pattern.  Figure  19,  shows  a  maximum 
gain  at  about  42  degrees  from  the  z  axis.  This  agrees  well 
with  the  38  degree  maximum  typically  observed  in  experiments 
[8:92]  . 

Figures  20  through  23  show  the  difference  patterns  of 
the  thin  wire  model  for  different  planes  of  observation. 
Instead  of  looking  in  the  X-Z  plane  (phi*0),  these  patterns 
are  taken  45  degrees  from  the  X-Z  plane  and  90  degrees  from 
the  X-Z  plane,  which  is  the  Y-Z  plane.  These  figures 
indicate  the  patterns  in  the  Y-Z  plane  (Fig.  22  and  23)  are 
almost  exactly  the  same  as  those  in  the  X-Z  plane  (Fig.  18 
and  19).  But  the  patterns  at  phi  »  45  degrees  are 
different.  The  pattern  at  nhi  ■  45  degrees  for  theta 
polarization  is  shaped  the  same  as  that  for  phi  >0  or  90 
degrees.  However,  the  theta  polarization  pattern  at  45 
degrees  is  shifted  about  6  dB  lower  than  the  other  patterns. 
This  suggests  that  a  three  dimensional  view  of  the  theta 
component  of  the  difference  pattern  is  a  "warped  donut"  with 
valleys  between  the  principal  planes.  Unfortunately,  this 
illustrates  the  antenna  has  less  than  perfect  circular 
symmetry . 

The  pattern  at  45  degrees  for  the  phi  polarization 
(Fig.  20)  is  about  2  dB  higher  than  its  counterparts.  This 
"reverse  shift"  also  indicates  there  is  a  relative 
polarization  breakdown  in  the  difference  pattern  as  the 
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Figure  23. 


Figures  24  through  27  show  the  sub  and  difference 
patterns  at  2  GHz  for  the  wire  grid  model.  The  maximum  gain 
of  the  sum  patterns  for  the  phi  and  theta  polarizations  are 
2.582  dB  and  3.545  dB,  respectively.  This  is  close  to  the  3 
dB  maximum  linear  polarization  gain  expected  from  a  spiral 
antenna.  The  beamwidths  (3  dB)  for  th  phi  and  theta  sum 
patterns  are  71.0  and  64.7  degrees,  respectively. 
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Figure  24 


r>R  1 D  1Q0EL  < 


0H2..  SUM  PATTERN 


Overall,  the  wire  grid  model  results  are  very  similar 
to  the  thin  wire  model  results.  There  is  a  slight 
descrepancy  in  the  phi  polarization  difference  patterns. 
The  wire  grid  result  (Fig.  26),  which  is  taken  in  the  X-Z 
plane,  shows  a  measurable  gain  on  the  z  axis.  This  behavior 
is  unexplained. 

Figures  28  through  31  and  32  through  35  show  more  wire 
grid  model  patterns  at  1  GHz  and  3.75  GHz,  respectively. 


RID  1S0F.L 


UM  PATTERN 


Figure  34 


Figure  35. 


Another  interesting  aspect  of  radiation  study,  beside 
the  radiation  patterns  themselves,  is  to  observe  the 
currents  developed  on  the  antenna.  In  this  study,  the 
currents  have  been  plotted  along  one  arm  of  the  thin  wire 
model.  The  actual  plots  show  the  expansion  mode  numbers 
versus  the  coefficients  of  the  current  expansion.  The 


expansion  mode  numbers  are  merely  an  indexing  of  the 
piecewise  sinusoidal  functions  that  go  from  1  to  99  as  the 


arm  is  traversed  from  the  feed  to  the  outer  tip.  To 
understand  the  meaning  of  the  plots,  the  mode  number  must  be 


converted  into  something  meaningful.  This  is  done  in  Table 
V.  This  table  converts  the  mode  number  into  the 
circumference  about  the  antenna  in  which  the  particular 
expansion  function  is  centered. 


Table  V 

Mode  Circumference  Conversion 


Mode  Number  Circumference  (cm) 


1 

34 

67 

2.50 

5.83 

13.53 

2 

3  5 

68 

2 . 57 

5.97 

13.88 

3 

36 

69 

2.63 

6.12 

14.24 

4 

37 

70 

2.70 

6.28 

14.61 

5 

38 

71 

2.77 

6.44 

14.99 

6 

39 

72 

2.84 

6.61 

15.38 

7 

40 

73 

2.92 

6.78 

15.78 

8 

41 

7. 

2.99 

6.96 

!6.  19 

9 
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10 
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3.23 
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17.48 

12 

45 

78 
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13 

46 
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14 

47 
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8.  12 

18.87 

1  5 

48 
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3.58 

8.32 

19.36 

IS 

49 

82 

3.67 

8.54 

19.86 

17 

50 

83 

3.77 

8.76 

20.38 

18 

51 

84 

3 . 86 

8.91 

20.91 

19 

52 

85 

3.96 

9.22 

21.48 

20 

53 

86 

4.07 

9.46 

22.01 

2  1 

54 

87 

4.17 

9.70 

22.57 

22 

55 
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4.28 

9.95 

2  3.16 

23 

56 
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4.39 

10.21 

23.76 

24 

57 
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4.50 

10.48 

24.38 

25 

58 

91 

4.62 

10.75 

25.01 

26 

59 

92 

4.74 

11.03 

25.65 

27 

60 

93 

4.86 

11.31 

26.32 

28 

61 

94 

4.99 

11.61 

27.00 

29 

62 

95 

5.  12 

11.91 

27.70 

30 

63 

96 

5.25 

12.22 

28.42 

31 

64 

97 

5.39 

12.53 

29. 16 

32 
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98 

5.53 

12.86 

29.91 

33 

99 

5.67 

13.  19 

30.69 
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Figures  36  through  40  show  the  amplitude  of  the  current 
coefficients  (normalized  to  one)  at  various  frequencies  for 
sum  pattern  generation.  These  currents  are  taken  from  the 
four  arm  thin  wire  model. 


SUM  M0De.  AT  S03  1HZ 


Figure  36. 


As  Mosko  [9]  indicated,  there  should  be  current 
resonance  where  the  circumference  of  the  antenna  is  one 
wavelength  for  sum  patterns.  Figure  38,  the  2  GHz  case, 
shows  resonance  around  circumferences  of  15  cm,  21  cm,  and 
29  cm.  The  wavelength  in  this  case  is  15  cm.  There  are  two 
extra  resonance  areas  on  the  spiral.  Extra  resonance  peaks 
also  occur  in  currents  forming  difference  patterns.  These 
currents 


t 


from  the  thin  wire  model,  are  plotted  in  Figures 


Mosko  [9]  explained  there  should  be  current  resonance 
where  the  circumference  equals  two  wavelengths  in  the 
difference  pattern  generation.  This  has  also  been  supported 
from  near  field  measurements  done  by  John  Dyson  at  the 
University  of  Illinois.  However,  Figure  43,  the  case  of  2 
GHz,  shows  resonances  at  circumferences  of  7  cm,  14  cm,  20 
cm,  and  28  cm.  The  only  resonance  should  occur  at  a 
circumference  fo  30  cm.  This  is  probably  where  the  last 
resonance  peak  "wants"  to  be,  but  it  must  shift  to  allow  the 
current  to  die  on  the  tip  of  the  arm. 

Figures  46  through  54  show  the  phases  of  current 
coefficients  producing  sum  and  difference  patterns  at 
various  frequencies.  They  show  that  the  phase  of  the 
current  changes  rapidly  in  low  amplitude  areas.  Again, 
these  plots  are  taken  from  the  currents  formed  in  the  four 
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Figure  54. 


The  currents  on  the  arms  of  the  wire  grid  model  cannot 
be  plotted  as  easily  as  those  in  the  thin  wire  model  because 
many  expansion  functions  overlap  transverse  and  longitudinal 
segments  of  wire.  However,  it  was  noted  in  this  model  that 
the  expansion  functions  that  contained  cross  hatches  had 
very  insignificant  coefficients,  except  near  the  center  of 
the  antenna.  This  would  help  explain  why  the  wire  grid 
difference  patterns  differed  more  with  the  thin  wire  results 


0 


than  the  sum  patterns.  This  is  because,  in  general,  the 
bulk  of  the  difference  pattern  current  is  closer  to  the 
center  of  the  antenna  than  the  bulk  of  the  sum  pattern 
current . 

The  extra  resonance  areas  in  the  current  distributions 
may  best  be  explained  by  an  inadequacy  in  the  modelling. 
The  open  circuited  tips  on  the  ends  of  the  arms  cause  a 
reflection  of  the  outward  travelling  current.  What  results 
is  a  standing  wave  pattern.  To  correct  this  situation,  it 
has  been  suggested  that  a  re  si  stive- induct ive  load  be  placed 
on  the  tips  of  the  arms  to  absorb  the  outward  travelling 
energy  [16].  To  do  this,  a  50  +  j50  ohm  load  was  placed  in 
the  thin  wire  model  on  the  point  connecting  the  outer  two 
segments  of  each  arm.  The  resulting  current  distributions 
for  the  sum  and  difference  pattern  generation  at  2  GHz  are 
given  in  Figures  55  and  56. 
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Figure-  55. 
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Figure  56 


There  is  very  little  change  in  the  current 
distributions  for  the  loaded  and  unloaded  thin  wire  models. 
Perhaps  a  much  larger  impedance  would  show  more  improvement. 
There  is  also  very  little  change  in  the  loaded  radiation 
patterns  from  the  patterns  shown  previously.  Figures  5? 
through  60  give  these  patterns  in  the  X-Z  plane  at  2  GHz. 
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Figure  58. 
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Figure  60. 


Results  for  radar  cross  section  are  not  given  in  this 


report 


Chapter  5 


Conclusions  and  Recommendations 

The  ESP  code  produced  radiation  patterns  very  close  to 
what  was  expected,  even  though  the  models  produced  a 
standing  wave  current  distribution  on  the  antenna  arms. 
This  standing  wave  pattern  results  from  the  models' 
inadequacy  in  absorbing  the  outward  travelling  current  on 
the  tips  of  the  arms. 

Since  there  are  so  few  differences  in  the  results  with 
the  thin  wire  grid  modelling,  it  is  only  necessary  to  use 
thin  wire  modelling  to  explore  radiation  patterns  or 
currents  of  interest.  However,  it  is  still  necessary  to  use 
patch  modelling,  which  is  not  feasible  on  a  VAX  11/750 
computer,  to  explore  the  structural  mode  characteristics  of 
radar  cross  section. 

Further  work  may  be  done  to  experiment  with  various 
loads  on  the  ends  of  the  arms  to  attempt  to  isolate  a  sigle 
area  of  current  amplitude  resonance  in  a  radiating  arm. 
This  load  should  be  the  complex  conjugate  of  the  impedance 
seen  by  the  outward  travelling  current  at  the  tip  of  each 
arm.  This  load  will  change  with  frequency,  as  the  impedance 
seen  at  the  tip  changes  with  frequency.  The  load  will  also 
affect  the  bandwidth.  In  particular,  it  may  raise  the  lower 
frequency  cutoff. 

If  a  single  load  cannot  eliminate  the  current 
reflection,  which  will  certainly  occur  when  a  large 
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frequency  range  is  tested,  another  solution  may  be  to 


attempt  a  series  of  loads  near  the  ends  of  each  arm.  This 
may  further  restrict  the  antenna  operating  bandwidth. 

Regardless  of  how  the  loading  is  attempted,  the 
elimination  of  current  reflection  would  probably  produce 
more  accurate  radiation  patterns.  This  is  because  actual 
spiral  antennas  usually  radiate  far  enough  from  the  outer 
circumference  so  that  the  tip  reflection  is  negligible.  The 
patterns  from  a  correctly  loaded  model  may  better 
demonstrate  over  what  angular  sector  the  frequency 
independence  holds.  Also,  the  radar  cross  section  may  be 
checked  to  see  if  loading  on  the  arm  tips  validate  the 


theory  given  in  Chapter  1 


Appendix  A 


Thin  Wire  Model  Input 

The  input  defining  the  thin  wire  model  geometry  for  the 
ESP  code  can  be  given  by  a  subroutine.  The  ESP  code  gives 
the  option  of  reading  the  rectarqular  coordinates  of  the 
wire  segment  endpoints,  or  generating  them  through  a 
subroutine  called  WGEOM.  The  option  to  use  subroutine  WGEOM 
is  determined  by  a  variable  known  as  IRGM  in  the  first  read 
statement  of  the  main  program.  If  this  variable  is  zero, 
the  program  uses  WGEOM  to  define  the  wire  geometry  [13:32]. 

The  goal  of  subroutine  WGEOM  is  to  define  the  following 
variables . 

IA ( I )  ■  an  array  (dimensioned  by  the  total  number  of 

segments)  which  defines  endpoint  A  of  the  Ith 
segment  with  a  point  number 

IB ( I )  »  an  array  (dimensioned  by  the  total  number  of 

segments)  which  defines  endpoint  B  of  the  Ith 
segment  with  a  point  number 

X , Y , Z  «  arrays  (dimensioned  by  the  total  number  of 

points)  which  define  the  rectangular 
coordinates  of  the  point  location  in  meters 

NM  ■  total  number  of  segments 

NP  *  total  number  of  points 

NAT  ■  total  number  of  wire  to  plate  attachments 

NWG  -  number  of  wire  geometries  -  this  variable  is  passed 
from  the  main  program  to  subroutine  WGEOM 

VG ( I )  *  complex  voltage  at  the  A  endpoint  of  the  Ith 

segment 

ZLD(I)  -  complex  impedance  at  the  A  endpoint  of  the  Ith 
segment 
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WV  *  wavelength  -  this  variable  is  passed  from  the  main 
program  to  subroutine  WGEOM 

In  addition,  the  array  variables  NS A ,  NPLA ,  VGA,  BDSK, 

and  ZLDA  exist  for  a  general  subroutine  WGEOM,  but  need  not 

be  defined  unless  wire  to  plate  attachments  are  involved  in 

the  geometry.  Also,  the  variables  NFS  1  and  NFS 2  exist  and 

need  to  be  defined  as  zero  when  mutual  coupling  between  feed 

points  is  neglected,  as  in  this  study. 

The  subroutine  used  to  define  this  input  information 

for  the  four  arm  thin  wire  model  is  listed  as  follows. 

SUBROUTINE  WGEOM  ( IA  ,  I B  ,  X  ,  Y  ,  Z  ,  NM  ,  NP ,  NAT  ,  iIS A  ,  NPLA  ,  VGA  , 

* BDSK, ZLDA, NWG , VG , ZLD , W V , NFS  1 ,NFS2 ) 

DIMENSION  IA< 1 )  , IB( 1  )  ,X{ 1  )  ,Y( 1  )  ,Z(  1  )  ,NSA( 1 ) ,NPLA(  1  )  , 
*BDSK( 1 ) 

COMPLEX  VGA (  1  )  , ZLDA (  1  ) , VG (  1  )  , ZLD ( i ) 

NM  *  404 
NP  *  405 

C  DEFINE  POINTS  FOR  ARM  #1 

DO  63  1-1,100 
Z(I)  »  0.0 
FF  -  REAL ( I } 

PHEE  -  (FF-1.0)*10. 

PHEEP  »  PHEE*3. 14159265/180. 

PH 90 0  «  PHEE/900  . 

R  «  0 . 3  9  8  *  (  10.**PH900)*0.01 
X(I)  -  R*COS ( PHEEP  > 

83  Y ( I )  *  R*SIN ( PHEEP ) 

C  DEFINE  ENDPOINTS  FOR  SEGMENTS  FOR  ARM#  1 

DO  84  1-1,99 
IA( I  )  «  I 

84  IB ( I  )  -  1  +  1 

C  DEFINE  POINTS  FOR  ARM #2 

DO  85  1-101,200 
Z { I )  «  0.0 
FF  -  REAL ( I ) 

PHEE  -  ( FF- 101. )*10. 

PHEEP  »  (PHEE+90. )  *  3 . 14159265/180. 

PH 90 0  -  PHEE/900. 

R  -  0 . 398* < 1 0 . ++PHS00 ) *0. 01 
X(I)  -  R*COS ( PHEEP ) 

85  Y ( I )  -  R*SIN ( PHEEP ) 

DEFINE  ENDPOINTS  FOR  SEGMENTS  FOR  ARM #2 
DO  86  1-100,198 
IA ( I )  »  1+1 


C 


o  o  o 


86  IB ( I )  =  1+2 

C  DEFINE  POINTS  FOR  ARM #3 

DO  87  1*201,300 
Z(I)  -  0.0 
FF  *  REAL ( I  ) 

PHEE  =  (FF-201  .  )*10. 

PHEEP  *  (PHEE+180. )*3. 14159265/180. 
PH900-PHEE/900. 

R  *  0 . 398* ( 10.**PH900)*0. 01 
X(I)  -  R*COS ( PHEEP ) 

87  Y { I )  *  R*SIN ( PHEEP ) 

C  DEFINE  ENDPOINTS  FOR  SEGMENTS  FOR  ARM #3 

DO  88  1-199,237 
I A ( I )  -  1+2 

88  IB ( I )  -  1+3 

C  DEFINE  POINTS  FOR  ARM#4 

DO  89  1-301,400 
Z(I)  -  0.0 
FF  -  REAL ( I ) 

PHEE  -  (FF-301  .  ) *10. 

PHEEP  -  iPHEE+270. )*3. 14159265/180. 

PH 90  0  -  PHEE/900. 

R  -  0. 398*( 10.**PH900 )*0. 01 
X(I)  -  R*COS ( PHEEP ) 

89  Y ( I )  »  R*S IN ( PHEEP ) 

C  DEFINE  ENDPOINTS  FOR  SEGMENTS  FOR  A RM # 4 

DO  90  1-298,396 
IA ( I )  -  1+3 

90  IB ( I )  -  1+4 

C  INDICATE  NO  ATTACHMENTS  OR  COUPLING 

NAT  -  0 
NFS  1  -  0 
NFS2  -  0 


DEFINE  FEED  STRUCTURE  FOR  SUM  PATTERN 


X  (  40  1  ) 

Y  (  4  0  1  ) 

Z  (  40  1  ) 

X  (  4  0  2  ) 

Y  (  4  0  2  ) 

Z  (  4  0  2  ) 

X ( 403  ) 

Y ( 4 03  ) 

Z  (  4  0  3  ) 

X ( 404  ) 
Y( 404  ) 

Z ( 404  ) 

X ( 40  5  ) 

Y ( 40 5  ) 

Z  ( 405  ) 
IA(  397  ) 
IL'(  397  ) 


0. 

0. 

0. 

0.00199 

0. 

0. 

0. 

0.00199 

0. 

-0.00199 

0. 

0. 

0. 

«  -0.00199 

»  0. 

-  1 
-  402 
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IA ( 398  ) 

as 

402 

IB( 398  ) 

S3 

401 

IA ( 399 ) 

s* 

101 

IB ( 399 ) 

S3 

403 

IA ( 400  ) 

SB 

403 

IB ( 400  ) 

m 

401 

IA ( 40  1  ) 

=* 

201 

IB  (  40  1  ) 

Si 

404 

IA ( 402  ) 

as 

404 

IB ( 402  ) 

SB 

401 

IA ( 40 1 ) 

m 

301 

IB ( 403 ) 

m 

405 

IA ( 404  ) 

m 

405 

IB ( 404 ) 

m 

401 

VG ( 398 ) 

m 

( 1 . 0, 0. 0) 

VG (400  ) 

m 

o 

• 

o 

• 

o 

> • 

VG ( 40  2  ) 

m 

1 

_A 

e 

O 

o 

« 

o 

**— * 

VG (404  ) 
RETURN 

END 

m 

o 

• 

o 

as 

1 

e 

O 

When  difference  pattern  feeding  is  desired,  the  last 
six  lines  of  the  subroutine  can  be  substituted  with  the 


following. 

VG ( 3  98  ) 
VG (400) 
VG  (402  ) 
VG ( 404  ) 
RETURN 
END 


«  (1.0, 0.0) 

-  (-1. 0,0.0) 
-  (1.0, 0.0) 

-  (-1.0, 0.0) 


When  a  complex  impedance  of  5Q+j50  ohms  is  placed 
between  the  last  two  segments  of  each  arm,  the  following 
statements  must  be  added  before  the  RETURN  statement. 


ZLD ( 99 ) 
ZLD ( 19A) 
ZLD ( 2  97  ) 
ZLD( 396  ) 


(50.0,50.0) 

(50.0,50.0) 

(50.0,50.0) 

(50.0,50.0) 


81 


n  o 


Appendix  B 


Wire  Grid  Model  Input 


A  subroutine  WGEOM  is  written  to  provide  the  input  data 

for  the  geometry  of  the  wire  grid  model.  An  explaination  of 

the  purpous  and  a  general  description  of  a  subroutine  WGEOM 

is  given  in  Appendix  A  and  need  not  be  repeated  here.  The 

following  is  a  listing  of  the  WGEOM  fortran  subroutine  for 

the  wire  grid  model,  with  comments  to  aid  in  understanding 

the  accomplishments  of  each  section  of  code. 

SUBROUTINE  WGEOM ( IA , IB , X , Y , Z , NM , NP , NAT , NSA , NPLA , VGA , 
*BDSK , ZLDA , NWG , VG , ZLD,WV,NFS1 ,NFS2  ) 

DIMENSION  IA(1),IB(1),X<1),Y(1),Z{1) ,NSA( 1 ) ,NPLA( 1 ) , 

*  BDSK ( 1  ) 

COMPLEX  VGA (1 ) , ZLDA ( 1 ) , VG ( 1 ) , ZLD ( 1 ) 

C  DEFINE  INNER  AND  OUTER  EDGE  POINTS  FOR  SPIRAL# 1  WITH 

45  DEGREE  INCREMENTS.  TUB  FIRST  23  POINTS  ARE  INNER. 
THE  SECOND  23  POINTS  ARE  OUTER. 

DO  150  1-1,23 
FF  -  REAL ( I ) 

PHEE  «  ( FF- 1 . ) *45. 

PHEEP  -  PHEE*3. 14159265/180. 

BB  -  0.001 08*PHEE 
RR  »  (  1t).**BB)*0.395*0.01 
X ( I )  -  RR*COS ( PHEEP ) 

Y ( I )  -  RR*SIN ( PHEEP ) 

Z ( I )  -  0. 

150  CONTINUE 
DO  151  1-24,46 
FF  -  REAL ( I ) 

HEE  -  (  FF-24.  )*45. 

PHEEP  »  PHEE*3. 14159265/180. 

BB  »  0.00106*PHEE 
RR  -  ( 10. **BB)*0. 4437*0. 01 
X(I)  ■  RR*OOS ( PHEEP ) 

Y ( I )  «  RR*SIN (PHEEP ) 

Z( I )  ■  0. 

151  CONTINUE 

C  SEGMENTATION  FOR  ARM#  1  •  THE  FIRST  22  SEGMENTS  FORM 

C  THE  INNER  EDGE.  SEG.  23-44  FORM  THE  OUTER  EDGE.  SEG. 

C  45-67  FORM  CROSS  HATCHES. 

DO  152  1=1,22 
IA ( I )  -I 
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IB ( I )  *  1+1 

152  CONTINUE 

DO  153  1=23,44 
IA ( I )  *  1+1 
IB ( I )  -  1+2 

153  CONTINUE 

DO  154  1*4  5,67 
IA ( I )  *  1-44 
IB ( I )  *  1-21 

154  CONTINUE 

DEFINE  INNER  AND  OUTER  EDGE  POINTS  FOR  SPIRAL#2. 
POINTS  47-69  ARE  INNER.  POINTS  70-92  ARE  OUTER. 

DO  155  1*4  V , 69 
FF  *  REAL ( I ) 

PHEE  *  { FF-47 . ) *45. 

PHEEP  »  (PHEE+90. )*3. 14159265/180. 

BB  -  0.001 08*PHEE 

RR  -  { 10. **BB)*0. 395*0. 01 

X(I)  -  RR*COS (PHEEP ) 

Y ( I )  -  RR*SIN ( PHEEP ) 

Z ( I )  -0. 

155  CONTINUE 

DO  156  1-70,92 
FF  «  REAL ( I ) 

PHEE  -  (FF-70. ) *4 5 . 

PHEEP  -  (PHEE+90 . )*3. 14159265/180. 

BB  *  0  .  0  0  1  0 6 *PH EE 

RR  -  ( 10. **BB)*0. 4437*0. 01 

X ( I )  -  RR*COS (PHEEP ) 

Y(I)  -  RR+SIN ( PHEEP ) 

Z ( I )  -  0. 

156  CONTINUE 

SEGMENTATION  FOR  ARM #2 .  SEG. 68-89  FORM  THE  INNER 
EDGE.  SEG. 90-111  FORM  THE  OUTER  EDGE.  SEG. 112-134 
CROSS  HATCHES. 

DO  157  1-68,89 
IA ( I )  *  1-21 
IB ( I )  -  1-20 

157  CONTINUE 

DO  158  1-90,  1 1  1 
IA ( I )  -  1-20 
IB ( I )  -  1-19 

158  CONTINUE 

DO  159  1*112,134 
IA ( I )  »  1-65 
IB ( I )  *  1-42 

159  CONTINUE 

DEFINE  INNER  AND  OUTER  EDGE  POINTS  FOR  SPIRAL#3. 
POINTS  93-115  ARE  INNER,  116-138  ARE  OUTER. 

DO  160  1-93,115 
FF  *  REAL ( I ) 

PHEE  «  (FF-93. ) *45. 

PHEEP  -  (PHEE+180. )*3. 14159265/180. 
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BB  -  0.001 08*PHEE 

RR  *  <  10.**BB)*0.395*0.01 

X  ( I  )  -  RR*COS  (  PHEEP  ) 

Y ( I )  *  RR*SIN (PHEEP ) 

Z(I)  -  0. 

160  CONTINUE 

DO  161  1-116,138 
FF  -  REAM  I) 

PHEE  »  (FF-1 16. ) *45. 

PHEEP  *  (PHEE+180 . ) *3. 14159265/180. 

BB  «  0.001 06*PHEE 

RR  »  (10. **BB)*0. 4437*0.  01 

X(I)  -  RR*COS (PHEEP ) 

Y ( I )  »  RR*S IN ( PHEEP  ) 

Z(I)  -  0. 

161  CONTINUE 

SEMENTATION  FOR  ARM#3.  SEG. 135-156  FORM  THE  INNER 
EDGE.  SEG. 157-178  FORM  THE  OUTER  EDGE.  SEG. 179-201 
ARE  CROSS  HATCHES. 

DO  162  1-135,156 
IA ( I )  -  1-42 
IB ( I )  -  1-41 

162  CONTINUE 

DO  163  1-157,  178 
I A ( I )  «  1-41 
IB ( I )  -  1-40 

163  CONTINUE 

DO  164  1-179,201 
IA ( I )  -  1-86 
IB ( I )  -  1-63 

164  CONTINUE 

DEFINE  INNER  AND  OUTER  EDGE  POINTS  FOR  SPIRALI4. 
POINTS  139-161  ARE  INNER,  162-184  ARE  OUTER. 

DO  165  1-139,161 
FF  -  REAL ( I  ) 

PHEE  »  ( FF- 139. ) *45. 

PHEEP  ■  (PHEE+270 .) *3. 14159265/180. 

BB  -  0.001 08*PHEE 

RR  -  ( 10.**BB)*0.395*0.01 

X(I)  -  RR*COS (PHEEP ) 

Y ( I )  »  RR*SIN (PHEEP ) 

Z(I)  -  0. 

165  CONTINUE 

DO  166  1-162, 184 
FF  «  REAL ( I ) 

PHEE  »  (FF-162. ) *4 5 . 

PHEEP  -  (PHEE+270. )*3. 14159265/180. 

BB  -  0.001 06*PHEE 

RR  *  ( 10. **BB)*0. 4437*0. 01 

X(I)  -  RR*COS ( PHEEP ) 

Y ( I )  »  RR*SIN (PHEEP) 

Z(  I )  -  0  . 

166  CONTINUE 


o  o 


C  SEGMENTATION  OFR  A  RM#4 .  SEG.  202-223  FORM  THE  INNER 

C  EDGE.  SEG. 224-245  FORM  THE  OUTER  EDGE.  SEG. 246-268 

C  ARE  CROSS  HATCHES. 

DO  167  1-202,223 
IA ( I )  -  1-63 
IB ( I )  -  1-62 

167  CONTINUE 

DO  168  1-224,245 
IA ( I )  -  1-62 
IB ( I )  -  1-61 

168  CONTINUE 

DO  169  1-246,268 
IA ( I )  -  1-107 
IB ( I )  ■  1-84 

169  CONTINUE 

THE  ARMS  AND  SEGMENTS  ARE  DEFINED. 

THE  REST  IS  FOR  THE  FEED  STRUCTURE. 

X ( 185)  -  0. 

Y( 185)  -  0. 

Z( 185)  -  0. 

X( 186)  -  0.0019 
“f 186)  -  0. 

Z(  186  )  -  0. 

X ( 1 87  )  -  0. 

Y( 187)  -  0.0019 
Z(  187)  -  0. 

X (  188  )  -  -0.00  19 
Y(  188)  -  0. 

Z ( 788)  -  0. 

X ( 1 89 )  -  0. 

Y ( 189)  -  -0.0019 
Z ( 189)  «  0. 

IA ( 269 )  -  1 
IB ( 269  )  -  186 
IA ( 270 )  -  186 
IB ( 270  )  -  185 
IA ( 2 7 1 )  -  47 
IB  (  2  7  1  )  -  187 
IA  (  2  72  )  -  187 
IB ( 2 72  )  -  185 
IA ( 273 )  -  93 
IB ( 2  73 )  -  188 
IA ( 274 )  «  188 
IB ( 2  74 )  -  185 
IA  (  2  75  )  *  139 
IB ( 2  75  )  -  189 
IA ( 2 76  )  -  189 
XB (  2 76  )  -  185 
NM  «  276 
NP  »  189 
NAT  «  0 
NFS  1  »  0 
NFS 2  «  0 
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o  o 


THE  FOLLOVJING  FOUR  DELTA  GAP  GENERATORS  ARE  FOR  SUM 
PATTERN  FEEDING. 

VG ( 270 )  -  (1.0, 0.0) 

VG  (212)  -  (0.0, 1.0) 

VG (274)  -  (-1.0, 0.0) 

VG  (276)  «  (0.0, -1.0) 

RETURN 
ENL 

When  difference  pattern  feeding  is  desired,  the  last 

six  lines  of  the  subroutine  can  be  su’»stituted  with  the 

following . 

VG (270)  »  ( 1.0, 0.0) 

VG (  272  )  -  (-1.0, 0.0  ) 

VG ( 2  74  )  «  ( 1 .0,0.0) 

VG ( 2  76  )  *  (-1.0, 0.0  ) 

RETURN 

END 
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Appendix  C 

Plate  Model  Programming 


In  a  program  involving  plates,  the  geometry  input  may 
be  defined  by  allowing  the  main  program  of  the  ESP  code  to 
read  the  rectangular  coordinates  of  each  corner  of  each 
plate.  However,  when  many  plates  are  used  in  a  model  and 
equations  can  be  used  to  generate  these  coordinates,  it  is 
simpler  to  use  a  subroutine  to  give  the  geometry.  Such  a 
subroutine  can  be  used  by  a  slight  alteration  of  the  main 
program  of  the  ESP  code.  The  following  is  a  section  of  the 
main  program  before  the  alteration. 

TOUCH  -  0.001*WV 

IOK  -  1 

READ(5,*)NPLTS 

IF ( NPLTS «  EQ. 0 )  GOTO  462 

IOKT  -  1 

DO  464  NPL-1, NPLTS 

READ { 5 , * )  NCNRS(NPL) , SEGM ( NPL ) , I REC ( NPL ) , IPN(NPL)  , 

*IGS (NPL ) 

DO  466  NCNP-1 ,NCNRS (NPL) 

READ ( 5 , * ) PCN ( 1 , NCNP , NPL ) , PCN ( 2 , NCNR , NPL ) , 

*PCN { 3 , NCNR, NPL ) 

466  CONTINUE 

CALL  PLPLCK { PCN,IPL,ICN,NCNRS(NPL) , TOUCH , NPL , IOK ) 

IF(IOK.EQ.0)IOKT-0 
464  CONTINUE 

IF( IOKT.EQ. 0 )G0T09374 

The  READ  statement  in  the  third  line  is  the  seventh 
READ  statement  in  the  main  program.  The  variable  being 
read,  NPLTS,  is  the  total  number  of  plates  in  the  geometry. 
If  the  number  of  plates  is  aero,  the  program  goes  to  line 
462,  which  occurs  just  before  the  tenth  READ  statement. 

The  loop  denoted  by  DO  464  is  used  to  read  the 
following  variables  [13:36]. 


NCNRS(NPL)  »  the  number  of  corners  in  plate  NPL 

SEGM(NPL)  3  maximum  segment  size  of  a  surface  patch  monopole 
(1/2  of  expansion  function)  in  wavelengths, 
typically  0.25 

IREC(NPL)  *  plate  indicator  for  plate  NPL;  it  is  zero  when 
plate  is  polygonal,  one  when  plate  is 
rectangular 

IPN(NPL)  ■  polarization  indicator,  chosen  to  be  three  when 
both  polarizations  of  expansion  functions  are  on 
the  plate 

IGS(NPL)  »  number  of  generating  side  in  a  later  subroutine 
for  mode  expansion;  it  is  chosen  to  be  zero  when 
the  longest  side  is  used 

The  loop  denoted  by  DO  466  is  used  to  define  the  X,  Y, 
and  Z  coordinates  for  each  corner  of  plate  NPL.  Then 
SUBROUTINE  PLPLCK  is  used  to  check  the  plate  geometry  to 
make  sure  the  corners  of  each  plate  are  in  the  same  plane 
[13:106].  Finally,  the  variables  I0K  and  XOKT  are  used  a 
flag  variables  which  are  used  to  abort  the  program  when 
SUBROUTINE  PLPLCK  sends  the  main  program  the  abort  message. 
The  GOTO  9374  sends  the  program  to  the  end  of  execution. 

To  modify  the  program  to  accept  a  plate  defining 
subroutine,  a  new  variable  must  b  read  by  the  main  program 
to  indicate  that  a  subroutine  called  WPLATE  is  being  used. 
This  new  variable  is  named  NPIND .  It  may  assume  the 
following  values. 

NPIND  -  0  for  plates  defined  conventionally  with  DO  464  and 
DO  466  loops 

NPIND  -  1  for  plates  defined  by  SUBROUTINE  WPLATE 

NPIND  must  be  set  to  zero  when  there  are  no  plates.  The 

same  passage  of  the  ESP  main  program  with  the  program 
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modifications  is  given  below.  New  and  changed  statements 
are  indicated  with  stars  in  the  left  hand  column. 

TOUCH  -  0.001*WV 

*  IOKT  -  1 
IOK  ■  1 

*  READ ( 5 , * )  NPLTS,  NPIND 
IF(NPLTS.EQ.O)  GOTO  462 

*  IF(NPIND.EQ.O)  GOTO  2891 

*  CALL  W PLATE  (NPLTS  , NCNRS ,  SEGN ,  IREC  ,  T  t'  -f ,  IGS ,  PCN  ) 

*  DO  3245  NPL-1, NPLTS 

*  CALL  PLPLCK(PCN,IPL,ICN,NCNRS(NPL)  TOUCH , NPL , IOK ) 

*  IF( IOK.EQ.O )IOKT-0 
*3245  CONTINUE 

*  GOTO  3333 
IOKT  «  1 

*2891  DO  464  NPL-1, NPLTS 

READ ( 5 ,  *  )  NCNRS(NPL) ,SEGM(NPL) ,)  ,J(NPL) , IPN (NPL) , 

* IGS ( NPL ) 

DO  466  NCNR-1 , NCNRS(NPL) 

READ (5,*)  PCN ( 1 , NCNR, NPL ) , PCN ( 2 • NCNR, NPL ) , 

*PCN(3, NCNR, NPL) 

466  CONTINUE 

CALL  PLPLCK ( PCN, I PL, ICN, NCNRS ( M V L ) , TOUCH, NPL, IOK) 

IF( IOK.EQ.O )IOKT-0 
464  CONTINUE 

*3333  IF(IOKT.BQ.O)  GOTO  9374 

At  this  point,  it  is  necessary  to  point  out  that  the 
corner  number  NCNR  should  be  placed  on  either  consecutively 
clockwise  or  counterclockwise  corners  around  a  plate.  If 
not,  SUBROUTINE  PLPLCK  will  reject  the  geometry.  Also  note 
that  the  only  variable  passed  from  the  main  program  to 
SUBROUTINE  If  PLATE  is  the  number  of  plates,  NPLTS. 

The  following  is  a  sample  SUBROUTINE  WPLATE  that 
generates  the  plates  for  one  arm  of  the  antenna  with  45 
degree  polar  angle  incrementation,  as  in  wire  grid 
modelling.  The  number  of  plates,  NPLTS,  is  22  to  model  the 
full  arm. 

SUBROUTINE  WPLATE (NPLTS , NCNRS , SEGM, IREC , IPN, IGS, PCN) 
DIMENSION  NCNRS ( 30 ) , SEGM( 30 ) , IREC ( 30 ) , IPN ( 30 ) , IGS ( 30 ) , 
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*PCN (3,4,30) 

C 

C  THE  DIMENSIONING  ALLOWS  FOR  A  MAXIMUM  OF  30  PLATES  AND 

C  FOUR  CORNERS  ON  EACH  PLATE. 

C 

DO  30  1*1, NPLTS 
NCNRS(I)  -  4 
SEGM(I)  »  0.25 
IREC(I)  *  0 
IPN(I)  *  3 
IGS(I)  *  0 
30  CONTINUE 

DO  40  1-1, NPLTS 
FF  -  REAL ( I ) 

PHEE  -  (FF-1 . ) *45 . 

PHEEP  -  PHEE*3. 14159265/180. 

PHEE0  -  PHEE+45. 

PH  -  PHEE0*3. 14159265/180. 

BB  -  0.001 08*PHEE 

RR  -  (  10.**BB) *0. 395*0. 01 

CC  -  0. 00108*PHEE0 

SS  -  ( 10.**CC)*0.395*0.01 

PCN (1,1,1)  -  RR*COS ( PHEEP ) 

PCN ( 2 , 1,1)  -  RR*SIN (PHEEP) 

PCN (3,  1,1)  -  0. 

PCN (  1,2,1)  -  SS*COS ( PH ) 

PCN (2,2,1)  -  SS *S IN ( PH  ) 

PCN (3, 2, I)  -  0. 

DD  -  0.00106* PH EE0 

TT  -  { 10. **DD)*0. 4437*0. 01 

EE  -  0. 00106* PH EE 

UU  -  ( 10. **EE)*0. 4437*0. 01 

PCN ( 1 , 3 , 1 )  -  TT*COS ( PH ) 

PCN (2,3,1)  -  TT *SIN ( PH ) 

PCN (3,3,1)  •  0. 

PCN (1,4,1)  «  UU*COS( PHEEP) 

PCN ( 2 , 4 , I )  -  UU*SIN ( PHEEP ) 

PCN (3,4,1)  *  0 . 

40  CONTINUE 
RETURN 
END 

This  subroutine  can  be  used  to  run  the  CPU  execution 
time  experiment  mentioned  in  Chapter  3.  Simply  change  the 
variable  NPLTS  read  in  the  seventh  READ  statement  of  the 
main  program. 

If  a  geometry  requires  both  wires  and  plates,  both 
SUBROUTINE  WPLATE  and  WGEOM  may  be  used.  The  wire  to  plate 
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attachments  would  be  put  into  SUBROUTINE  WGEOM 


C1 
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